Effect of bond lifetime on the dynamics of a short-range attractive colloidal system 
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We perform molecular dynamics simulations of short-range attractive colloid particles modeled 
by a narrow (3% of the hard sphere diameter) square well potential of unit depth. We compare the 
dynamics of systems with the same thermodynamics but different bond lifetimes, by adding to the 
square well potential a thin barrier at the edge of the attractive well. For permanent bonds, the 
relaxation time r diverges as the packing fraction <f> approaches a threshold related to percolation, 
while for short-lived bonds, the ^-dependence of r is more typical of a glassy system. At intermediate 
bond lifetimes, the ^-dependence of r is driven by percolation at low cf>, but then crosses over to 
glassy behavior at higher <j>. We also study the wavevector dependence of the percolation dynamics. 

PACS numbers: 82.70.Gg, 61.20.Lc, 64.60.Ak, 82.70.Dd 



I. INTRODUCTION 

Colloidal systems, in which particles are dispersed in a 
fluid, have an enormous relevance in industrial applica- 
tions, owing to the possibility of chemically or physically 
tuning the interaction between the particles and the re- 
sulting possibility of designing materials with novel prop- 
erties From the point of view of basic research, 
colloidal systems are playing a very important role in the 
development of the physics of liquids, since they open 
up significantly the range of values of physically rele- 
vant parameters. For example, novel phenomena arise 
when the range of particle-particle interaction becomes 
significantly smaller than the size of the particle or when 
the system is composed of colloidal particles with signif- 
icantly different size or mobility. 

An interesting phenomenon that is often observed in 
colloidal suspensions, but is absent from atomic or molec- 
ular systems, is particle clustering and gelation. The gel 
is an arrested state of matter at small values of the pack- 
ing fraction, a non-ergodic state capable of supporting 
weak stresses. The formal description of gel formation 
in colloidal systems has been receiving considerable at- 
tention recently 0, IE IE S 13 • Recent numerical work 
has also focused on the gelation problem [E El • Inter- 
esting studies have attempted to provide formal connec- 
tions between the formation of a gel and the formation of 
a glass, both being disordered arrested states of matter. 
It is not a coincidence that such theoretical studies fo- 
cus on colloidal systems, where colloid-colloid interaction 
can be finely tuned, allowing for a careful test of theo- 
retical predictions. Indeed, colloids appear to be ideal 
systems for unraveling the physics of gel formation. Un- 
derstanding the key features of the interaction potential 
that stabilize the gel phase will probably have an impact 
also on our understanding of the protein crystallization 



problem where the possibility of generating crys- 

tal structures is hampered by the formation of arrested 
states. 

Sterically stabilized colloidal particles provide an ex- 
perimental realization of a system in which the particle- 
particle interaction can be well modeled by the hard 
sphere potential [TjJ. When this is the case, addition 
of many small non- adsorbing polymers leads, due to de- 
pletion mechanisms, to an effective short-range attrac- 
tion between the colloidal spheres |T1 . Il5| . Neglecting 
the effects of the solvent on the dynamics of the colloidal 
particles, and integrating out the behavior of the smaller 
polymers, one has an experimental realization of a short- 
range potential, with a tunable short-range attraction 
between particles. The size of the small polymers tunes 
the range, while their concentration controls the strength 
of the attraction uq. 

At high packing fraction cf> w 0.6, these colloidal sys- 
tems exhibit the usual hard sphere glassy dynamics. 
When the range of interaction is smaller than about 10 
percent of the hard-sphere diameten the glas s transition 
line can show re-entrant behavior |1. Htj IitI HH flij. I20L 
I2H l22l 123 . |24| . That is, in a particular range of cj), the 
liquid can be brought to structural arrest by either in- 
creasing or decreasing the ratio T/uq, where T is the 
temperature. Experimentally, dynamical arrest phenom- 
ena in short-range attractive colloids are observed not 
only at high density, as discussed above, but also in the 
low packing fraction region. In this case, the arrested 
material is commonly named a gel. The gel state dis- 
plays peculiar features like the appearance of a peak in 
the static structure factor, for very large length scales (of 
the order of several particle diameters), that is stable in 
time, as well as a non-ergodic behavior in the density cor- 
relation functions and a finite shear modulus [2{| . These 
solid-like, disordered, arrested features have prompted 
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the appealing conjecture that these colloidal gels can be 
viewed as the low-density expression of the high-density 
glass line, with both phenomena being driven by the same 
underlying mechanism of arrest 0- E3- HE- However, 
such a connection between gelation and the attractive 
glass is non-trivial, as pointed out in Ref. [2jj . The pres- 
ence of an intense pre-peak in the static structure factor 
has suggested also the possibility that, in colloidal sys- 
tems, the gel phase is related to a phase separation pro- 
cess |28|, l2E HE HH H3, H3 ■ Indeed, hard sphere systems 
with short-range attraction added tend to phase separate 
into a colloid rich phase (liquid) and a colloid poor (gas) 
phase. Whether the interaction between this phase sepa- 
ration and the re-entrant glass line can bring about a gel 
phase via arrested phase separation [ssl ] is an idea which 
is also under current investigation. At very low T, diffu- 
sion limited cluster aggregation [3 HE HE H3 may be 
another way to irreversibly obtain a clustered gel (and a 
frozen pre-peak in the structure factor). 

In the present study, as a step in the process of un- 
derstanding gelation in colloidal systems in the absence 
of phase separation, we focus on the dependence of 
the dynamics on a purely kinetic factor, the lifetime of 
the particle-particle bond. We introduce a Hamiltonian 
model of a short-range attractive colloid, for which we 
can tune the bond lifetime, without affecting the ther- 
modynamics. We have been inspired by the recent work 
of Del Gado and coworkers |3q|. where a lattice model 
was introduced to study the influence of bond lifetime 
on the slow dynamics of gelling systems. Here, we model 
a colloidal system as an ensemble of particles interact- 
ing with a short-range square well, a model sufficiently 
realistic to properly describe the physics of short-range 
systems, but at the same time ideal for studying particle 
bonding and percolation since a bond is unambiguously 
defined by the limits of the square well. In particular, 
we study the interplay between percolation and the glass 
transition and find that there is a crossover from a perco- 
lation dominated regime, to one controlled by the glass 
transition. Differing from Ref. pjsl ]. we find that above 
the percolation threshold, changing the lifetime of bonds 
merely rescales long time behavior of the system, leaving 
intact the glassy a-relaxation. Furthermore, we explore 
the long lifetime case with regards to the dependence of 
observing percolation on the wavevector used to probe 
the system, making contact with experimental observ- 
ables. Finally, as a contribution towards clarifying the 
differences between gels and glasses, we study the same 
model in the limit of permanent bonds, where percolation 
becomes the relevant arrest process in the system. 



II. MOLECULAR DYNAMICS SIMULATIONS - 
THE MODEL 

We perform collision-driven molecular dynamics sim- 
ulations of a binary mixture of particles interacting 
through a narrow square well pair potential. Although 



colloidal systems are more properly modelled using Brow- 
nian dynamics, we use molecular dynamics due to its ef- 
ficiency in the case of step-wise potentials. While the 
short-time dynamics is strongly affected by the choice 
of the microscopic dynamics, the long term structural 
phenomena, in particular close to dynamical arrest, are 
rather insensitive to the microscopic dynamics [39j |. We 
use a 50:50 binary mixture of 700 particles of mass m 
with diameters <jaa = 1.2 and ubb — 1 (setting the unit 
of length). The hard core repulsion for the AB interac- 
tion occurs at a distance gab = {&AA + c\b.b)/2. The 
depth of the well uq is 1, and the width Ay of the square 
well attraction is such that Ay /(cry + Ay) = 0.03 for all 
interactions between particles of type i and j. T is mea- 
sured in units of uq, time t in (JBBim/uo) 1 ^ 2 ■ This system 
has been extensively studied previously [2E HE EE Ell • 

The phase diagram of this system is reproduced from 
Ref. [33j in Fig. ^ F° r this model, both the dynamical 
arrest curves and the spinodal curve have been calcu- 
lated. The glass line (determined by extrapolating the 
diffusion coefficient calculated in simulation to zero ac- 
cording to a power-law 14111 shows both a repulsive and 
an attractive glass branch j^E ■ The numerical glass lines 
are well described by the predictions of mode-coupling 
theory (MCT) |42l | , after an appropriate mapping is per- 
formed |4lll43| . Fig. Q] also reports the static percolation 
line (defined as the locus of points in (<fi, T) such that 50% 
of the configurations possess a spanning, or percolating, 
cluster of bonded particles) and the estimated location 
of the liquid-gas spinodal (the locus of T below which 
spinodal decomposition occurs in simulation). 

It is important to note that in this model the attractive 
glass line ends on the spinodal line on the large 4> branch, 
proving that arrested states at low (j) in this model can 
arise only as a result of interrupted phase separation |33| . 
It also confirms that, if the MCT predictions for the lo- 
cation of the attractive glass are not properly rescaled in 
the (j> — T plane, an incorrect location of the glass line 
with respect to the spinodal line is predicted. 

In order to study the effect of bond lifetime on the dy- 
namics, wc add to the edge of the square well an infinites- 
imal barrier of tunable height h (see Fig. [2J , thereby sta- 
bilizing bonds formed when particles become trapped in 
the attractive square well of the pair potential [44J . As 
the barrier is infinitesimal, the portion of phase space oc- 
cupied is negligible, and hence the thermodynamics of the 
system is unaffected. For numerical reasons, in the code 
we have implemented a barrier width of 3 x 10~ 4 cjbb 
checking that the static structure of the system is not 
affected by this tiny but non-zero width. 

To the extent that the thermodynamics are unaffected 
by the barrier, configurations drawn from equilibrium 
simulations of the h — case are also representative con- 
figurations of the system when h ^ 0. Thus, results for 
different h > are obtained using 30 or more indepen- 
dent initial configurations obtained by equilibrating the 
system for h = 0. This technique alleviates the compu- 
tational burden when working with large values of h. 
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FIG. 1: Phase diagram of the square well binary mixture (re- 
produced from [33]), showing the percolation line (squares), 
approximate location of the spinodal (dashed line), repulsive 
and attractive glass transition lines (solid line). State points 
studied here are shown as filled circles. The highlighted state 
points (T = 0.5, cj> = 0.52 and T = 1.5, <f> = 0.595) refer to 
those presented in Fig. [(J and Fig. |7| respectively. 



V(r) 



u f 



FIG. 2: Schematic of the pair potential. Here shown are the 
hard sphere core diameter a, narrow square well of depth 
Mo = 1 and width A, with A/(A + a) = 0.03. The height 
h of the barrier controls the bond lifetime, and hence the 
microscopic dynamics of the system. 



In this article we focus on the time dependence of the 
(collective) density-density correlation function (dynamic 
structure factor). The dynamic structure factor, the cor- 
relation function typically accessed in scattering exper- 
iments, is given by F q (t) = (p q (t)p- q (0)) /S(q), where 
p i^> = ^E^i ex P(-^-^), S(q) = (M0)|) 2 is the 
static structure factor, (.) denotes an ensemble average, 
fi is the position vector of a particle, q is a wavevector and 
i labels the N particles of the system. We also make use 



of the correlation function for type A particles only, de- 
fined similarly as F A (t) = (p q (t)p A q (0)) /S A (q), where 



i=l ex P (- 



iq 



r t ), s A ( q ) = (\(4(o)\y 



the partial static structure factor for type A particles, 
and the summations are over the Na particles of type A. 
The qualitative behavior of F q {t) and F A (t) is the same 
for this system. 



III. THE INFINITE BOND LIFETIME CASE 

We start discussing the case of bonds of infinite life- 
time, i.e., the case where h — > 00. In this limit, a well de- 
fined model for continuum percolation is generated. The 
spatial distribution of the particles is fully determined by 
the equilibrium properties of the square well potential 
(and hence static correlations are known and precisely 
defined) while the dynamics is the dynamics of a sys- 
tem constrained by irreversible bonds. The possibility of 
generating equilibrium structures with h = to be used 
as starting configurations for the case h 7^ allows us 
to completely decouple issues arising from the bond life- 
time from issues associated with non-equilibrium prop- 
erties and aging also when the packing fraction is larger 
than the percolation value. The averaging over different 
starting configurations allows us to properly sample con- 
figuration space. To study the h — > 00 case we perform 
simulations at T/uq — 0.5, for about 30 different values 
of packing fraction, as indicated in Fig. ^ Numerically, 
we achieve the infinite limit by setting h = 1000, a value 
high enough so that we never observe bond breaking. 

The 4> dependence of P(4>), the fraction of particles 
belonging to the spanning cluster, provides a way of de- 
tecting the location of the percolation point. In all per- 
colated configurations, we observe the presence of only 
one spanning cluster. When finite size effects are negligi - 
ble, P ~ \(j)-(j) p \ where (3 is a critical exponent 45, 46] . 
Fig. |21 shows P{<j>). The arrow in the figure indicates 
(j) p = 0.23, which we identify as the value of the packing 
fraction at which a spanning cluster is found in 50 per- 
cent of the configurations [47j ■ To estimate the effect of 
bonding on the dynamics, we show in Fig.^Jthe packing 
fraction and wavevector dependence of F q (t). 

For <p < <f) p (Fig. EJa)), correlation functions decay 
to zero, independently of the value of q, as expected 
for a system where only diffusive clusters of finite size 
are present. For (f> > <f> p (Fig. ^b)), an "infinite" span- 
ning cluster is present. On increasing <f> the size (mass) 
of the spanning cluster increases progressively, incorpo- 
rating 90% of the particles in the system already when 



(f> = 0.32 (see Fig. |2Jl. For <js > <f) p , wavevectors are char- 
acterized by correlation functions which do not decay to 
zero any longer, reflecting the presence of a non-relaxing 
component of the density fluctuations. 

Close to percolation, only for very small wavevectors 
does F q (t) go to a non-zero plateau of height f q , also 
called the non-ergodicity factor. On increasing packing 
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FIG. 3: P as a function of </>■ The arrow indicates the (f) p = 
0.23 at which 50% of the sampled configurations contain a 
percolating cluster. 



fraction, the amplitude of the plateau increases signifi- 
cantly, as shown in Fig.^Jc). Simultaneously, correlation 
functions at larger and larger wavevectors show a finite 
non-zero long time limit. The wavevector and <f> depen- 
dence of f g is shown in Fig. |3] 

The appearance of a non-zero f q , whose amplitude and 
width grow on increasing <j> is consistent with the onset of 
a percolation transition and the loss of ergodicity of the 
particles in the infinite spanning cluster. The inverse of 
the half- width of f q provides an estimate of the associated 
localization length. On increasing (f> beyond percolation, 
such a length decreases from infinity (or from the simu- 
lation box length in a finite size system) down to the di- 
mension of the particles, in analogy with the progressive 
decrease of the connectness length of the spanning clus- 
ter 0, ^3 . Owing to the large localization length close 
to percolation and to the tenuous structure of the per- 
colating cluster, close to percolation a non-zero f q value 
can be clearly detected only at very small wavevectors. 
On increasing <f>, the increase in the number of particles 
in the infinite cluster and the associated decrease of the 
length generate an increase in the amplitude and width 
of f q , making possible the numerical observation of a 
non-zero f q even at large q. 

It is interesting to compare the behavior of the non- 
ergodicity factor observed in the case of percolation with 
the case of the glass. The most striking difference is in 
the change of f q across the glass and percolation tran- 
sitions. In the case of glasses, f q shows a discontinuous 
jump, while in the case of percolation it increases from 
zero continuously. In the language of MCT , the per- 
colation transition is analogous to what is called a "type 
A" transition while the ordinary glass transition is of 
"type B" . 



Another important aspect is the fact that the width 
in q space of f q (e.g., the q value at which a curve in 
Fig. |S{a) reaches half its height) is of the order of the 
inverse of the nearest neighbor distance in the case of 
glasses (or even larger in the case of attractive glasses) 
while it is extremely small close to percolation. Only 
when most of the particles are part of the spanning clus- 
ter, does the width of f q become compatible with the 
one of glasses. This change in the width of f q reflects the 
significant difference in localization length of the parti- 
cles (the length scale on which particles are trapped in 
chiefly vibrational motion) at the glass transition (of the 
order of the nearest neighbor or of the bond distance) 
as opposed to the very large localization length at the 
percolation transition where a tenuous (almost massless 
in the thermodynamic limit) spanning cluster appears. 

It is important to note that in the case of glasses, the 
glass transition is marked by the arrest of density fluctu- 
ations on every length larger than the nearest neighbor 
distance, while in the case of percolation, the observation 
of a non-ergodic transition is strongly dependent on the 
observation length. To clarify this point, Fig.[3Jb) shows 
the <j) dependence of f q at several q values (i.e. a cut of 
the data shown in Fig.|S^a) at fixed q). We note that the 
steep increase of f q from zero occurs at larger and larger 
4> values on increasing q. This suggests that experiments 
— capable of measuring a non-zero f q with a finite pre- 
cision — restricted to a fixed q value will notice a loss of 
ergodicity in the sample, as reflected by a non-zero long 
time limit of the correlation function, only at a 4> value 
which may be significantly larger than the percolation 
packing fraction. 

In analogy to the connection between </> p and the van- 
ishing of P(<p) (also shown in Fig|Sfb)), one may define 
an apparent <j) c {q) based on the packing fraction at which 
the f q curves shown in Fig. [3b) cross a fixed value, con- 
trolled by the precision of the experimental technique 
in detecting a non-zero f q value. This <p c (q) could be 
considered an indicator of the percolation transition as 
observed at a particular q vector. 



IV. FINITE BOND LIFETIME: EFFECT OF 
BARRIER HEIGHT ON F q (t) 

When h has a finite value, the lifetime of the bond is fi- 
nite. Hence, a new time scale enters into the description 
of dynamics in the model. In particular, we are inter- 
ested in the modification of the density correlation func- 
tions introduced by the finite bond lifetime, and in the 
competition between the bond time scale and the caging 
time scale close to the glass transition. In Fig. a) we 
plot F g A (i) at <f> = 0.52 and T = 0.5 for qa B B ~ 2tt, for 
several increasing values of h. When h = 0, the decay 
of the correlation function does not show signatures of 
two-step relaxation, owing to the location of the state 
point in the re-entrant liquid portion of the phase dia- 
gram (state point is highlighted in Fig. 20]. As h 
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FIG. 4: Dependence of F q (t) on <f) and g for h — oo. (a) 
F q (t) at <p = 0.225 (just below the percolation threshold of 
0.23). Here, q<JBB = nqo for various n shown in the legend, 
with go = ?r/£ = 0.2408 (L is the length of the simulation 
box). At n — 2, it is more difficult to reduce noise in the 
data because of the small number of g-vectors available for 
averaging, (b) F q (t) at <f> = 0.38 (well above percolation), 
with go = 0.2867. (c) F q (t) at qo b b ~ 27r for various <f>. 
From such curves we determine the plateau height f q . 
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FIG. 5: (a) Plateau height f q as a function of g for various 
<j>, for infinite h (irreversible permanent bonds). Legend gives 
<j> values, (b) Connection to percolation: f q as a function 
of <j> for various q. The legend indicates values of qcr b b for 
f q curves corresponding to filled symbols. The curve labeled 
with open circles shows P, the average fraction of particles 
participating in a percolating cluster (taken from Fig. [JJ . 



increases, two new features appear: (i) a slowing down of 
the correlation function and (ii) the emergence of a two 
step-relaxation process. Correlation functions decay to a 
rather high plateau before decaying to zero at long times. 
When h = oo, the correlation function does not decay to 
zero any longer. In Fig. |SJb), despite the emergence of a 
plateau, the long time behavior is merely rescaled with 
respect to the h — case. The shape of the correlation 
functions does not change significantly in the long time 
regime and indeed curves for different (finite) h values 
can be superimposed on the h = curve, as shown in 
Fig-Elk). As a suitable scaling time we choose the time 
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FIG. 6: (a) The dynamic structure factor at qoss = 6.37 
for the state point cj> = 0.52, T = 0.5 (see Fig. for various 
values of the barrier height h. The attractive plateau persists 
longer as h increases, i.e., as the bond lifetime is increased. 
Panel (b) shows that the functions collapse onto a common 
curve at long times, when t is rescaled (see text). 



at which the correlation functions reach 20% of their ini- 
tial value. More precisely, the density correlation func- 
tions are plotted as a function of a rescaled time, t/th, 
where t h = t 20 (/i)/t 20 (0), and FffaQi)) = 0.20. 

We next discuss the case in which, in the absence of a 
barrier, the liquid is close to a repulsive glass transition, 
marked by the presence of a two-step relaxation in F^(t). 
We show in Fig. da) the correlation functions for differ- 
ent barrier heights for the state point T = 1.5, cf) — 0.595 
(highlighted in Fig.QJ. At this higher T, the h = sys- 
tem behaves as a hard sphere binary mixture. The slow 
dynamics is thus characteristic of repulsive glass dynam- 
ics and shows a well defined plateau. On increasing h, 
one observes a progressive modification of the correlation 
function at short times, and the emergence of the "gel" 
plateau (highlighted in the inset). When time becomes 




FIG. 7: (a) The dynamic structure factor at q<JBB = 6.67 for 
the state point cj> = 0.595, T = 1.5 (see Fig. for various 
values of the barrier height h. (b) We again see the stabiliza- 
tion of the attractive plateau (highlighted in the inset), and 
simple time rescaling at long times. 



comparable to the bond lifetime, the correlation func- 
tion leaves the "gel" plateau and approaches the caging 
plateau, following the same dynamics as in the h = 
case, as clearly indicated by the superposition of curves 
for different h values on a rescaling of the time (Fig.|Z£b)). 
This superposition indicates that the lifetime of the bond 
indeed acts to renormalize the "microscopic time". An 
increase in h increases the time required for breaking 
the particle-particle bonds which in turn increases the 
timescale of the breaking and reforming of cages. 

The fact that the slow relaxational processes are unaf- 
fected by the bond lifetime (apart from a trivial scaling 
factor depending on h) is particularly reassuring for theo- 
ries of the glass transition which connect static properties 
to dynamic properties, like MCT. The peculiarity of this 
model is indeed the fact that static properties are inde- 
pendent of h. Hence, according to MCT all dynamical 
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FIG. 8: The dependence of T20 on h and <j>, obtained from 
F^(t) at qctBB = 2-7T. The curves are for h = (circles), 
h = 1 (squares), h = 2 (diamonds), /i = 3 (up-triangles) , 
h = 4 (left-triangles) and permanent bonds (down-triangles). 
For permanent bonds, T20 diverges near <f> = 0.43. For smaller 
values of h, T20 at smaller <j> tracks the divergence at <j> = 0.43, 
but then crosses over to glassy dynamics with a divergence at 
higher <j>. 

properties associated with caging should be independent 
of h. The scaling observed in Fig.EJb) and Fig.[7Jb) sup- 
ports such a hypothesis. We also note that the results re- 
ported here differ from those reported in Ref. [3^ , where 
clustering induced by the bonds was considered to be sig- 
nificantly connected to the glass transition phenomenon. 
One possible explanation of such a difference may lie in 
the fact that in the study of Ref. (3^ , at odds with the 
present model, the bond lifetime is strongly coupled to 
the structure of the system. 



V. CROSSOVER FROM PERCOLATION TO 
GLASSY DYNAMICS 

We now focus on the <f> dependence of the characteristic 
time for different values of h. To quantify the character- 
istic time, we use r 2 o, the time at which F^{t) decays to 
a value of 0.2. In Fig. [S] we show the dependence of T20 
on h and </> for wavevector modulus q = 2w /obb- 

We see from Fig. [H] (and from Fig. ^) that T20 for the 
h = system appears to diverge near <f> > 0.6. For 
the case of permanent bonds, h = 00, T20 instead shows 
an apparent divergence at <p c ~ 0.43. As discussed in 
the previous sections, this divergence is a manifestation 
of the percolation transition which has taken place at 
<t> p = 0.23 (see Fig. EJb), near <j> w 0.43 for qasB = 27r). 
At intermediate values of h, the <f> dependence of T20 is 
highly non-trivial. A crossover from the percolation be- 
havior (h — 00) to the glass behavior (h = 0) takes place 
at a typical time controlled by the value of h. This is 
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FIG. 9: The dependence of T20 on h for various packing frac- 
tions, at qoBB = 2-7T. Above cf> c , the characteristic time tends 
to an Arrhenius behavior with h, i.e., T20 ~ exp (h/T). 

most clearly seen for h — 3 and h = 4 in Fig. |SJ Indeed, 
deviations from the h — 00 case are expected when the 
lifetime of the bond becomes shorter than T2o(h = 00). 
In another way, on time scales shorter than the bond life- 
time, the bonds appear to be permanent and the system 
behaves like the h — 00 case. 

To estimate the role of h in slowing down the dynam- 
ics we report T20 as a function of h for various isochores 
in Fig. We find that above (p c , T20 approaches an Ar- 
rhenius behavior with respect to h as h increases, i.e., 
^20(0,^) ~ <?(<^) exp (h/T) where g(4>) is function only 
of the state point chosen. This factorization of time 
allows us to clarify that the main effect of bonding is 
to redefine the microscopic time scale of the dynamics. 
The bond lifetime does not affect the properties of the 
(a-relaxation) slow dynamics on approaching the glass 
transition. 

It is interesting to state the connection of the present 
findings with an earlier work |4(jj|. where the barrier was 
used to extend the high-plateau to timescales associated 
with a-relaxation of the native system, in the packing 
fraction region where a glass-glass transition was ex- 
pected. The present data clearly show that the barrier 
does not affect (except for a rescaling of the microscopic 
time) the true a-relaxation dynamics. On the other 
hand, large barrier values bring into the window of exper- 
imental observation the intra-well dynamics and move to 
inaccessible regions the a-relaxation. Under these condi- 
tions, the decay of the correlation function in the experi- 
mentally accessible window is limited to the high-plateau 
f q value, which coincides with the attractive glass. 

This suggests the possibility of simultaneously study- 
ing gelation and glassy dynamics within the same system, 
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by focusing on different time scales. For example, in a 
system with transient bonds, observation on time scales 
much shorter than the bond lifetime would reflect the 
percolation dynamics, and hence gelation, while observa- 
tion on time scales much longer than the bond lifetime 
would yield results driven by glassy dynamics. 

VI. CONCLUSIONS 

In this paper we have introduced a simple model for 
studying continuum percolation in a system with well de- 
fined spatial correlations between the particles. In this 
model, bonding ambiguities are suppressed by the square 
well shape of the potential. While studying this model 
with permanent bonds, we have focused on the behav- 
ior of the density-density correlation function across the 
percolation transition, defined as the packing fraction 
at which a spanning cluster appears. We have found 
that the behavior of the density fluctuations is signif- 
icantly different from the one characteristic of super- 
cooled liquids and glasses. In the percolation problem, 
the non-ergodicity parameter increases continuously from 
zero and the localization length is extremely large and 
becomes comparable to the one observed in glasses (of 
the order of about 0.1 of the nearest neighbor or less) 
only well inside percolation. In this respect, percolation 
(with infinite lifetime bonds) and the glass transition are 
two distinct phenomena with distinct experimentally de- 
tectable signatures. We have also shown that in the case 
of percolation, since the range of wavevectors where non- 
ergodic behavior is observed grows with q on increasing 
(j) (for <j) > (ftp), experiments at fixed wavevector (due 
to their intrinsic finite resoluation) detect a non-ergodic 
transition at a packing fraction larger than <p p . In the 
case of glasses, the observation of an crgodic to non- 
ergodic behavior is essentially identical at all q values 
(except at very large q, describing self intra-cage mo- 



tion [13). 

The model allows us also to study the effect of the fi- 
nite bond lifetime while altering neither the structure nor 
the thermodynamics of the system. A comparison at dif- 
ferent bond lifetimes is thus performed on configurations 
which are characterized by the same particle-particle cor- 
relation. This study confirms the results recently re- 
ported by Del Gado and coworkers |38( for a lattice model 
concerning the existence of a crossover in the dynamical 
properties from a percolation controlled dynamics to a 
glassy dynamics on increasing cj>, when the lifetime of the 
bond is longer than the microscopic particle dynamics 
in the absence of a barrier. However, our results differ 
from Ref. (3^ in that we find that the bond lifetime acts 
essentially as a redefinition of the microscopic time and 
does not alter any feature of the slow dynamics and of 
the scaling laws approaching the glass transition. Still, 
the dynamics at times shorter than the a-relaxation time 
is strongly affected by the finite lifetime of the bond. 
The addition of the barrier, which increases the bond 
lifetime, extends the duration of the plateau character- 
istic of short-range attractive glasses (Fig. |7{a) inset). 
Since here we are in the liquid regime, the duration of 
the high-plateau is controlled by the (tunable) bond life- 
time. For times longer than exp(/i/T) the correlation 
function leaves this plateau and approaches the (lower) 
plateau associated with caging dynamics (Fig. Ufa)), and 
then finally relaxes to zero leaving the intrinsic slow 
long time dynamics of the system intact. 
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